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ABSTRACT 

Ch 

Aims, 51 Oph is one of the few young Be stars displaying a strong CO overtone emission at 2.3 microns in addition to the near 
infrared excess commonly observed in this type of stars. In this paper we first aim to locate the CO bandheads emitting region. Then, 
we compare its position with respect to the region emitting the near infrared continuum. 
£1*) Methods. We have observed 51 Oph with AMBER in low spectral resolution (R=35), and in medium spectral resolution (R=1500) 

centered on the CO bandheads. 

Results. The medium resolution AMBER observations clearly resolve the CO bandheads. Both the CO bandheads and continuum 
i-C emissions are spatially resolved by the interferometer. Using simple analytical ring models to interpret the measured visibilities, we 

find that the CO bandheads emission region is compact, located at 0. 15°qq 4 AU from the star, and that the adjacent continuum is 
q coming from a region further away (0.25°qq 3 AU). These results confirm the commonly invoked scenario in which the CO bandheads 

5^ originate in a dust free hot gaseous disk. Furthermore, the continuum emitting region is closer to the star than the dust sublimation 

radius (by at least a factor two) and we suggest that hot gas inside the dust sublimation radius significantly contributes to the observed 
jrt 2 yum continuum emission. 
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^ 1. Introduction 

CO 

ON The environment of young stars is formed by a dusty and 
gaseous circumstellar disk in which the first stages of planet for- 

\Q mation are thought to take place. Hence, building a comprehen- 

<^> sive picture of the physical mechanisms governing this forma- 

OO tion requires the understanding of the distribution and the evolu- 
tion of both the dust and gas species. The content and dynamics 
of the gaseous component, which dominates the total mass of 

•rH the disk, control the accretion-ejection processes and will de- 
fine the final architecture of forming planetary systems. On the 
other hand, the dust, responsible for the continuum infrared ex- 
cess, plays an essential role in determining the structure of the 
protoplanetary disk itself, since its opacity dominates in the con- 
tinuum. Dust grains are also the first building blocks that will 
grow into kilometric planetesimals and then into eventual plan- 
ets. 

To study the planet formation taking place in the very inner re- 
gion of the disk, one needs to combine milli-arcsecond high 
angular resolution techniques with high spectral resolution in 
order to directly probe the physical processes at stake both 
in the continuum and in the emission lines. In this respect, 
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near infrared spectro-interferometry that associates long base- 
line near infrared interferometry with spectroscopy appears per- 
fectly suited. This technique indeed offers a unique way to spa- 
tially resolve both the continuum infrared excess and emission 
line components, that is to locate the regions of emission inde- 
pendently of any prior assumption regarding the physical mech- 
anisms at play. In that sense, infrared spectro-interferometry has 
already provided important results about the origin of the hydro- 
gen line Bry. Different authors have thus located the Bry emis- 
sion region, showing that it can arise either from the accreting 
columns of gas falling onto the star ( |Eisner| 2007| > or from out- 
flowing winds in which ionized matter is ejected from the star 
dMalbet et an|2Tj07l|TItuili et al.|[2007bl|Kraus et al.|[2008l ). 
In this paper we are interested in exploring lines which are more 
direct probes of gas in the disk itself. We present the first inter- 
ferometric observations of the 2.3pm CO overtone emission in 
the (B9) Be star 5 1 Oph. 5 1 Oph is one of the very few young 
stars where this emission is strong enough to be observed with 



infrared spectro-interferometry ( |Thi et al.]|2005| [Berthoud et al. 



2007 ). This star exhibits also a strong Bry emission, which, as- 
suming a magnetospheric origin (Muzerolle et al. 2004), indi- 



cates an accretion rate of M = 1 - 2.10~ r M Q /yr ( |Garcia Lopez 
et al. 2006[ |Brittain et aL||2007[). The dust component is rather 



compact (Lagage & Pantin ( 1994); |Liu et al.|(|2007|l) and shows 
a strong silicate feature in emission (Bouwman et al. 2001|> but 
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Fig. 1. AMBER calibrated spectrum of 51 Oph around the 2-0 
and 3-1 bands of the CO overtone at 2.3microns. For comparison 
purposes is plotted (dashed line) the same spectrum measured 
with the TNG spectrograph (L. Testi, private communication). 
Note that we did not plot the TNG spectrum between the two 
bandheads because of irrelevant instrumental artifact. 



no PAH (Mee us~et al.| 2001|l, consisten t with a fiat and geomet 



rically thin dusty disk (Malfait et al. 1998). All the previous 



analysis are pointing towards a disk seen nearly edge-on. 



2. Observations and data reduction 

51 Oph (K = 4.3) has been observed with AMBER, first in 
medium spectral resolution (MR-K, R=1500) with the 8m Units 
Telescopes (UTs) of the VLTI in September 2006, then in low 
resolution (LR-K, R=35) with the 1.8m Auxiliary Telescopes 
(ATs) in September 2007. The log of the observations is pre- 
sented in Table [T] 

Spectrally dispersed observations: The observations in medium 
spectral resolution have been performed in the wavelength range 
[2280, 2340]^m around the CO overtone bandheads emission. 
On the AMBER detector are recorded both the dispersed pho- 
tometry from which the spectrum of our target is derived, and 
the dispersed interference pattern from which the visibilities as 
a function of the wavelength are computed. A reference source 
is required to calibrate both observables (i.e. the spectrum and 
the visibilities). 

51 Oph calibrated spectrum: In order to correct for the wave- 
length dependent telluric and instrumental features in the 
AMBER spectra, we used the observations of the calibration star 
HD 170499. This star has a K4III spectral type and has intrinsic 
photospheric absorption features that need to be removed be- 
fore using its observations to derive the response of the atmo- 
sphere and instrument. In particular, cool stellar photospheres 
show strong CO absorption longward of 2.9 pm (e.g. [Wallace 
& Hinkle"] ( |1997| l). To correct for the photospheric emission we 
constructed a template photospheric spectrum averaging three 
of the best quality spectra of stars with spectral type as close as 
possible to our own from | Wallace & Hinkle] ( |1997) : HR 1457 
(K5III), HR 3275 (K4.5III), and HR 6705 (K5III). The average 
spectrum (normalized to the continuum) was then multiplied by 
the spectral slope of a K4 III star using the Pickles ([T998) mod- 
els available on the ISAAC web pages at ESO to obtain the tem- 
plate. The template spectrum was then smoothed to the AMBER 
spectral resolution and used to remove the photospheric signa- 
tures in the calibrator spectrum and derive the response of the 
system. This response was then used to correct the spectrum of 



51 Oph. In Fig.[Tjwe show the final AMBER spectrum compared 
to a K band spectrum of 5 1 Oph obtained at the TNG telescope 
with the NICS NIR spectrograph (Isella, priv. comm.). 
51 Oph dispersed calibrated visibilities: We computed visibil- 
ities using the standard data reduction algorithms described in 
itulli et al. ( 2007a|l, previously applying the AMDC sofware 
( |Li Causi et al.| |2007[ ) to take away the spurious fringes of the 
detector that can bias the estimation of the visibility. A dras- 
tic fringe signal-to-noise ratio (SNR) selection using the 20% 
best frames was performed (for both the source and the calibra- 
tor) to insure reliable absolute calibrated visibilities (Tat ulli et] 
al. 2007a). These latter were obtained by dividing the 51 Oph 



raw visibilities by the reference ones, assuming a diameter of 



1 .28mas for the calibrator, according to the catalog of Merand et 



al. ( 2005 1. Finally, we binned the spectrally dispersed visibilities 



in the wavelength direction in order to increase the SNR. The re- 
sults are presented in Fig. [2] and the values of the visibilities in 
the continuum (i.e. average on the continuum spectral range) are 
given in Table [T] 

51 Oph low resolution visibilities: LR-K data were reduced fol- 
lowing the same procedure as the MR-K ones and then averaged 
spectrally to derive one single visibility of the continuum per 
baseline, whose values are plotted in Fig.[3]and given in Table[TJ 
as well. 



3. Modelling 

Our MR-K data resolve spectrally the two CO bands (although 
not the rotational components within each band). For the longest 
baselines UT2-UT4 and UT1-UT4, the emission is spatially re- 
solved at all wavelengths, though with rather high visibilities 
of 0.8. Our LR-K data are consistent with the MR-K ones, dis- 
playing high visibilities between 0.8 and 1 in the continuum at 
2.2pm. 

To derive more information, we need to compare the observed 
visibilities to model predictions. We adopt simple geometrical 
models and describe both the continuum and bandheads emit- 
ting regions as coming from rings of constant surface brightness. 
At each wavelength we assume that the emission is the sum of 
the unresolved stellar flux F s , of a continuum excess emission 
F c and of emission in the CO bandheads F/, with ratios that we 
derive from observations. In particular, the ratio F C /{F S + F c ) is 
derived from the spectral energy distribution (SED) photometric 
1988} , fitting a Kurucz spectrum ( |Kurucz 



data (Waters et al. 



[T9791 ) for the star with T t = 10000K, R t = 5.3R Q . We find 
F C I(F S + F c ) = 25% and FJ(F S + F c ) = 30% at 2.2pm and 
2.3pm, respectively. The ratio Fi/(F S + F c + Fi) is given by our 
normalized AMBER spectrum (see Fig.[T] left) in the CO band- 
heads and is about 18% at the peak of the 2-0 band and 19% at 
the peak of the 3-1 band. 

We compute two families of models, one with narrow ^jrings and 
one with broad rings extending from the stellar radius (in the 
following, uniform brightness). The parameters of the model are 
the inclination (/), the position angle (PA) of the system and the 
size of the bandheads (/?/) and continuum (R c ) emitting region, 
respectively. The visibility as a function of the wavelength is 
given by: 



F s + F c V c( R c U PA) + FfVf(R h i, PA) 



V A = 



F$ + F A C + Ff 



(1) 



1 the width of the ring is fixed at the typical value of AR/R = 20%, 
the results being poorly dependent on its value as long as AR <K R. 
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Date 


Telescopes 


Baseline (m) 


PAD 


Wavelength (jim) 


Resolution 


Calibrator 


Absolute visibility 


10-09-2006 


U1-U2-U4 


55/82/121 


34/91/69 


2.3 


1500 


HD 170499 


1.0+0.1/0.8±0.05/0.8±0.03 


05-09-2007 


G1-D0-H0 


69/60/71 


326/80/17 


2.2 


35 


HD163955 


0.8±0.1/0.85±0.1/1.0±0.1 


05-09-2007 


G1-D0-H0 


65/46/70 


345/89/24 


2.2 


35 


HD163955 


1.0±0.15/1.0±0.1/0.87±0.1 


06-09-2007 


G1-D0-H0 


64/31/67 


0/98/27 


2.2 


35 


HD 172051 


0.85±0. 12/0.95±0. l/0.95±0. 1 
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Fig. 2. MR-K calibrated visibilities as a function of the wave- Fig. 3. LR-K calibrated visibilities as a function of the projected 
lengths for the three baselines and overplotted best model (black baselines, for the three observations, and overplotted best model 
thick line), for the ring case. (lines). In each panel, colors red, blue,and green correspond to 

baselines Gl -DO, D0-H0, and Gl -DO respectively. 



Table 2. Best-fit model parameters. 

Model gT TQ PA (°) K (AU) R,(MJ) 

Uniform Brightness 08 85^ 129^ 0.43^ 0.36^$ 

Narrow Ring 0.8 82+* 5 126+j 5 0.24+g : gg 0.15°f 04 



4. Discussion 

4.1. Sizes of emitting regions vs. modelling 

The uniform brightness and the narrow ring are the two extreme 
cases of a unique model of uniform distribution in which the 
width of the emitting region is varying. We emphasize that this 
geometrical modelling defines the metric that allows us to un- 
equivocally link the measured visibilities to the spatial bright- 
ness distribution of these regions, through the Zernicke-Van 
Cittert theorem (Goodman 1985). It does not introduce assump- 
tions regarding the mechanisms at the origin of these emissions. 
We present in Table [2] the results for the two scenarios. All the 
intermediate models would also work, with intermediate values 
for the sizes. Note that the uniform brightness returns the largest 
size because emission is added at small spatial scales down to 
the stellar radius. The uniform brightness has no proper phys- 
ical meaning since neither the dust nor the CO gas can sur- 
vive so close to the star. However, it represents the character- 
istic size of the emission region. As this quantity is commonly 
used in infrared interferometry, it is given here for this purpose. 
Conversely the narrow ring model is sensitive to the location 
where most of the emission is arising from, usually the inner rim 
of this emission. In the following, we will analyze our results in 
the light of the narrow ring only. We emphasize that the size of 
the CO emitting region is smaller than the one of the continuum, 
independently of the choice of our model. 

4.2. Inclination and position angle of the system 

We find a position angle of 126^ 5 degrees and an inclination of 
i = 82^j 5 degrees. This inclination is in agreement with most 
of the previous studies. |Berthoud et al.| ( |2007[ ) have, however, 



where V$ is the visibility of the region emitting the continuum 
and Vf that of the region emitting the bandheads. 
Fitting together the continuum and bandheads models, we have 
performed a least-square minimization test on a large grid of pa- 
rameters varying both radii from to 3AU, the inclination from 
to 90° and the position angle from to 1 80° . The \ 2 function 
presents a clear and unambiguous minimum, for two reasons: (i) 
the high value of the visibilities imposing that we are in the first 
lobe of the visibility function, hence that the emission regions 
are rather compact and, (ii) the broad range of baseline orienta- 
tions which put strong constraints on the inclination and the po- 
sition angle of the system. The parameters of the best-fit model 
are summarized in Table [2] and the result of the fit is shown in 
Fig.|2]and[3] for the narrow ring model case. 
We note that both the bandheads and continuum components are 
unresolved for the smallest UT baseline shown in the top panel 
of Fig. [2] This suggests that the small variations of the visibil- 
ity detected in the CO bandheads are probably remaining biases 
of the data reduction, given the low SNR of the interferograms. 
However, considering the error bars of the measurements, this 
(possible) bias has only a very marginal influence on the sizes 
we derive and does not impact on their interpretation. 
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proposed an alternative scenario in which the object would have 
an inclination lower than 36°. Our measurements are positively 
ruling out this latter option. 

4.3. Locating the CO bandheads emission region 

For the first time, we are spatially resolving the CO bandheads 
emission region in 51 Oph. We find that most of the CO band- 
heads are emitted at a radius of R/ = 0.15_PLAU from the 
central star. This is fully supporting the scenario in which the 
CO bandheads emission is arising from the inner region of a 
hot gaseous disk in Keplerian rotation, between 0.15AU and 
0.35 - 0.53AU, as derived from spectroscopic modelling of the 
CO bandheads ( |Thi et al.[ [2003) |Berthoud etlE] |2007]) , and 
where the temperature is hot enough (> 2000K) to excite the 
first overtone bands at 2.3/im. More precisely, these authors, 
though they are fitting the CO bandheads with a broad ring, 
are pointing out that most of the emission comes from a re- 
gion close to the inner radius of their gaseous disk model at 
~0.15AU. This is in excellent agreement with the value we de- 
rive in the present analysis. We note that one single ring enables 
us to reproduce the visibility in the two bandheads, indicating 
that both the 2-0 and 3-1 bands come from the same physical re- 
gion. We also note that the CO emission region is inside the near 
infrared continuum one, the latter being located at a radius of 
R c = 0.24+^AU. This result corroborates that the CO is emit- 
ted from a dust free region, regardless of the true origin of the 
continuum (see Sect. |4.4| i. At such a distance from the star, the 
survival of the CO gas from photo-dissociation must however 
be considered. Thi et al. ( 2005 \ suggested that two mechanisms 
can circumvent the photo-dissociation process: the self-shielding 
of the CO molecules that oc curs if the column density is high 



van Dishoeck & Black 



(1988)), as 



enough (N(CO) > 10 15 cnr 2 , 
well as the C + OH — > CO + H chemical reaction, which insures 
a continuous replenishment of this molecule. 



4.4. The origin of the continuum near-infrared excess 

The continuum excess of 5 1 Oph is weak and, according to the 
shape of its SED (Waters et al. 1988), appears to come from 
a small and tenuous disk of dust which is optically thin at most 
wavelengths, as already inferred by Malfait et al. ( 1998 (. We per- 
formed SED fitting with the MCFOST dusty disk code of |Pinte 
|et al.| P006| ), and found the same results with a low amount of 
dust present in the disk (Md us t ~ 10~ 8 M G ). The result of our 
modelling is shown in Fig. [4] This model however requires a 
dust inner edge located at the dust sublimation radius, which in 
the case of 51 Oph would be at Rdust ~ 1-2AU for a typical 
grain size distribution of spherical particle^] and a dust subli- 
mation temperature of T^ ust ~ 1500/T. Obviously this model, 
with such a large inner rim, is not compatible with our mea- 
sured visibilities as shown by Fig. [5] and for which we calculate 
a. x 2 °f Studying the extreme case where only big grains 
are present would shrink the dust sublimation radius down to 
Rdust ~ 0.56AU, still at least two times further out than the dis- 
tance suggested by our data. Interestingly, a similar problem was 
evidenced by Leinert et al. ( |2004[ ) on the same star, the 0.5AU 
mid-infrared size of 5 1 Oph derived from their MIDI measure- 
ments being too small to be compatible with the shape of the 
SED. 

Several scenarios may however circumvent this apparent dis- 




X (urn) 

Fig. 4. Spectral energy distribution of 5 1 Oph taken from liter- 
ature (jWaters et al. 1988 ) (crosses) and superimposed best fit 



model (solid line) using MCFOST dusty disk model. The main 
parameters that are fitted are the inner radius Rdust = 1 .2 AU, the 
outer radius R out = 400AU scale height at 1AU H = 0.04, the 
flaring exponent /3 = 1.1, the grain size distribution (see text), the 
surface density exponent p — —1.25 and the tot al mass of dust 
Mdust = 10~ 8 M o . We refer to 
description of the modelling. 



Pinte et al. 



(2006) for a thorough 




MR-K data 
(continuum) 



LR-K data 



dn(a) oc a 3 7 da, ranging from a„ 



with optical constants from Mathis & WhifTen ( 1989 ) 



0.03 jim to a n 



1 mm, 



50 100 
Projected baseline (m) 

Fig. 5. MR-K continuum visibilities (top panel) and LR-K vis- 
ibilities as a function of the baseline and and overplotted 
MCFOST model arising from SED fitting (lines). The color code 
is the same than in fig.|2]and[3] We can see that this model, with 
inner radius equal to the dust sublimation radius, is incompatible 
with the interferometric data, with a resulting ;f 2 of ~10. 



agreement of the location of the dust evaporation radius with 
respect to our derived inner edge position. First, the gas inside 
the dust sublimation may help the dust to survive closer in by 
absorbing a fraction of the stellar radiation available to heat the 
dust. According to Muze rolle et al.| ( p004| ), the gas accreting onto 
the star becomes optically thick to stellar radiation for accre- 
tion rates higher than M ~ 10~ 7 M Q /yr. The accretion rate of 
1 - 2.10~ 7 M G /yr derived from Bry luminosity would actually 
place the gas in 5 1 Oph at the frontier between the optically thin 
and optically thick regimes. 

Another possibility comes from the rapidly rotating nature of 51 
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Oph. Given its measured rotational velocity vsin/ = 270km/s 



(Dunkin et al. 



19971), and knowing that this star is seen almost 



edge-on, 5 1 Oph would rotate at 90% of its critical velocity given 
by Vera = ^2GMJ(3R t ) 300km/s. Therefore, 51 Oph is likely 
elongated with a drop of its gravity (g e ff) from pole to equator, 
and subsequently of its effective temperature ( von Zeip"el) 1924) . 
The star being cooler at the equator than at the poles, the dust dis- 
tributed in the equatorial plane will be heated less efficiently and 
the evaporation radius moved closer to the star. From the grav- 
ity darkening law T e ff x geff (solid body approximation), the 
effective temperature of 5 1 Oph at the equator can be estimated 
to be of the order of 75% of that of the pole, hence lowering the 
sublimation radius by a factor of 2. 

Finally, one last explanation would be that the infrared excess is 
not entirely originating from dust thermal emission. The gas in- 
side the dust sublimation radius could substantially contributes 
to the near infrared energy balance. [Muzerolle et al.| p004 ) have 
shown that inner optically thick gaseous zones are indeed ex- 
pected to emit a large continuum excess in the near infrared 
through free-free emission. Interestingly, this effect has been re- 



cently shown to take place in Herbig Ae/Be stars ( |Eisner et al. 
[2007j|Isella eFaTl[2058l|Tannirkulam et al.||2008p . 



4.5. 51 Oph: a classical Be star? 

5 1 Oph appears to be a peculiar source in an unusual transitional 
state. In the frame of Herbig Ae/Be stars, its SED presents a 
near infrared excess which is too small to account for a classi- 
cal puffed-up inner rim. In the frame of /3-Pic like stars, 51 Oph 
is also lacking the far infrared-excess bump associated with the 
presence of an outer dusty disk, the inner disk being emptied by 
a potential forming planet (Malfait et al. 1998 1. Instead these 
authors suggest that 51 Oph is undergoing an alternative evolu- 
tion scheme, without forming planets. Furthermore, the presence 
of strong CO overtone emission bandheads also makes 5 1 Oph 
quite a puzzling case. Such emission requires large column den- 
sities of warm gas in order to produce detectable emission. Such 
large column densities are rare except in sources with the largest 
accretion rates (Najita et al. |2007 1. Interestingly, the CO over- 
tone has also been detected in the B9 star HD58647 (Berthoud 
2008 ), for which the SED profile is similar to that of 51 Oph 
( jMalfait et aT] [T998J. In HD58647 how ever, there is less CO 
emission but its infrared excess is s tronger. |B erthoud] ( |2008[ l thus 
concluded that both stars are most likely classical Be stars sur- 
rounded by massive gaseous disk, though seen at different evo- 
lutionary stages. This scenario is in agreement with their high 
rotational velocities (vsin / = 270km/s and 118km/s for 51 Oph 
and HD58647 respectively) and gives credit to the hypothesis 
that most of the near infrared continuum emission - if not all 
- is arising from the circumstellar gas. HD58647 has been ob- 
served with the Keck interferometer (Mon nler et al.| |2005| ) and 
these authors have derived a rather compact size of 0.4AU for 
its near infrared continuum emission, slightly larger as that of 
51 Oph for an equivalent luminosity (L = 260L Q for 51 Oph vs. 
L = 25OL for HD58647). This is again below - or at the very 
lower limit of - the dust sublimation radius, pointing towards 
the same gaseous origin for the continuum than that of 5 1 Oph. 
Hence, the very inner environment of this type of stars seems 
to follow an intriguing evolution scheme where the dust pro- 
gressively dissipates leaving behind a massive gas-rich, strongly 
accreting disk. As the hot dust is vanishing, the relative contri- 
bution of the gas to the continuum infrared excess increases, its 
region of emission moving closer to the star. How the circum- 



stellar dust disappears remains rather unclear. Repeated studies 
combining spectroscopic detection of emission bandheads (CO, 
hydrogen) and interferometric measurements in a large sample 
of stars would certainly help improve our understanding on how 
these disks are evolving and dissipating. 



5. Summary 

The principal results are summarized here: 

- using the AMBER/VLTI interferometer, we have spatially 
and spectrally resolved the CO bandheads emission at 2.3jum 

- the inclination of the system is constrained by our observa- 
tions, and we confirm that the object is seen nearly edge-on. 

- assuming simple ring models for the distribution of both the 
near infrared continuum and CO emissions, we locate the 
bandheads emission region at a distance of 0.1 5° ° y 4 AU from 
the star, inside the continuum emission region, reinforcing 
the scenario in which the CO bandheads originate from the 
inner part of a dust free gaseous disk. 

- the continuum emission is located at a distance of 
0.24^q q3 AU, at least two times closer to the star than the dust 
sublimation radius. This suggests that the inner disk of gas is 
playing an important role, whether by shielding the star light 
and thus allowing the dust to survive closer to the star, and/or 
by contributing to a non negligible part of the near infrared 
continuum emission. 
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